Has DAMA Detected Self-Interacting Dark Matter? 
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We consider a model in which dark matter consists of a large self-interacting component (SIDM 
particles) and a small component with strong interactions with baryons (SIMPs). We show that the 
DAMA annual modulation signal can be caused by SIDM particles interacting with SIMPs trapped 
in iodine nuclei in the DAMA/Nal detector. 
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The nature of dark matter (DM) remains a mystery, 
despite many years of efforts to understand it. It has 
been shown 0] that coUisionless dark matter correctly 
predicts the large scale structure of the universe. It is less 
clear, however, if coUisionless dark matter is consistent 
with observations on the galactic scale and smaller. Some 
studies have suggested that there are discrepancies with 
observations. It has been found that: 

• More satellite galaxies are predicted in simulations 
than are observed 0- 

• Simulations of disk galaxies contain less mass and 
angular momentum than is observed Q- 

• Simulated dark matter halos of galaxies have 
steeper density profiles compared to observations 
i. 

Dark matter consisting of particles with strong self- 
interactions (SIDM particles] have been suggested as a 
solution to these problems p, 1^ . The cross section be- 
tween SIDM particles, crad, needed is in the range 0,0^ 
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Here msioM is the mass of the SIDM particle. While 
some other studies have casted doubt on the self- 
interacting dark matter hypothesis and have suggested 
upper limits on add of the order of ~ 10~^^ cm^ 
™siDM/GeV, see e.g. 0, in a recent new study it is 
argued that Udd could be much higher, about 3.6 x 10~^^ 
cm^ msiDM/GeV. 

Cross sections between dark matter particles could 
thus be of the same order as for particles interacting 
with each other via the strong force. This has led some 
to speculate that DM particles could interact with both 
themselves and with baryons through the strong force 
1^. DM particles interacting with baryons with such 
strong interactions had been postulated earlier by Dover, 
Gaisser, Steigman and Wolfram and have been named 
SIMPs (Strongly Interacting Massive Particles). SIMPs 
have been proposed as a solution to the ultra high en- 
ergy cosmic ray problem [Toj and as an explanation for 



the absence of cooling flow s IT]J . Constraints on SIMP 
dark matter are given in Il2l Il3| . 

It has been shown [l4. ITsI that SIMPs can bind 
to heavy nuclei, forming anomalously heavy isotopes. 
Searches for such isotopes have been performed in gold 
and iron recently 16] and have come out negatively. Lim- 
its on the abundance of anomalous gold isotopes range 
from 6 X 10"^^ for light isotopes ~ 3 GeV heavier than 
a gold nucleus to 1 x 10~^ for heavy isotopes up to 1480 
GeV heavier than a gold nucleus. In 13] these limits 
were used to constrain the contribution of SIMPs to the 
total energy density of the Universe as a function of the 
SIMP mass. 

Currently, there are about 20 dark matter direct de- 
tection experiments underway. In these experiments one 
tries to detect recoils of nuclei caused by collisions with 
DM particles. Most of these experiments use sophis- 
ticated background discrimination techniques to try to 
extract a DM signal. An alternative way to extract a 
DM signal from the background is to measure the annual 
modulation in the dark matter flux caused by the Earth's 
rotation around the Sun [H 113. The DAMA/Nal ex- 
periment 0] has used this technique and is the only one 
to have reported a positive signal. The results of this ex- 
periment are controversial because other more sensitive 
searches have not detected nuclear recoils due to DM in- 
teractions. For WIMPs interacting with nuclei through 
spin independent interactions, the available parameter 
space is very limited The case of spin dependent 

interactions is studied in Inelastic dark matter 
and mirror matter j23l ] have also been proposed as an 
explanation for the DAMA results. 

In this article we consider a scenario in which SIMPs 
make out a small fraction of the dark matter density of 
the Universe. We assume that the remaining part of the 
dark matter consists of SIDM particles with almost no in- 
teractions with nuclei but with strong interactions, sim- 
ilar in magnitude to their self-interactions, with SIMPs. 
We shall show that dark matter direct detection experi- 
ments can detect signals from a SIDM particle colliding 
with a SIMP trapped in a nucleus and that this effect 
may be responsible for the annual modulation signal ob- 
served in the DAMA experiment. Since in this model 
the self-interactions of the SIDM particles don't matter. 
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we shall call them WIMPs and refer to this model as the 
WIMP-SIMP model. 

The DAMA/Nal set up consists of 9 scintillating thal- 
lium doped Nal crystals of 9.7 kg each. Electrons and nu- 
clei recoiling after a collision cause emissions of photons 
that are detected using photomultiplier tubes. Nuclear 
recoils give rise to a weaker photon signal than electron 
recoils. The ratio of the two signal strengths is called the 
quenching factor and is 0.3 for sodium recoils and 0.09 
for iodine recoils. The measured signals are expressed in 
electron equivalent recoil energies, i.e. the energy of a re- 
coiling electron needed to account for the observed signal. 
After a total exposure of 107731 kg. day an annual modu- 
lation in the energy range 2-6 KeV was observed with an 
amplitude of 0.0192 ± 0.0031 counts/day/kg/KeV 

In a WIMP-SIMP scenario, assuming that SIMPs 
don't bind to sodium, all recoils contributing to the an- 
nual modulation should come from anomalous iodine iso- 
topes. When a WIMP scatters off a nucleus, the nucleus 
recoils with an energy, ii^iccoih of: 
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where 9 is the scattering angle in the center of mass 
frame, v the velocity of the WIMP and /iwiMP, nucleus 
is the reduced mass of the WIMP-nuclcus system. The 
DAMA data have been analyzed assuming different types 
of interactions of nuclei with WIMPs and a wide variety 
of WIMP halo models For our purpose, the case of 
spin independent interactions is the most interesting. In 
this case the WIMP-nucleus cross section scales with the 
square of the atomic mass number and with the square of 
the relative WIMP-nucleus mass. This means that most 
collisions are with iodine nuclei. Instead of exploring the 
full window allowed by the DAMA data, we'll use that 
a WIMP mass of 50 GeV is consistent with the DAMA 
data for spin independent interactions. For A/wimp ~ 50 
GeV the factor (a* wimp, nucleus)^ /Afnucieus in © is ap- 
proximately 11 GeV. If we assume that the quenching 
factor for anomalous iodine isotopes is the same as for 
ordinary iodine, then constraining the WIMP and SIMP 
masses such that: 
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where /j.wimp,i+simp is the reduced mass of the WIMP- 
anomalous iodine system and Afj is the mass of an iodine 
nucleus, will guarantee that the WIMP-SIMP model is 
consistent with the DAMA data. 

We have calculated the WIMP-SIMP cross section 
in terms of the abundance of anomalous iodine iso- 
topes by assuming an isothermal halo model |l7j for 
WIMPs and the value of 11 GeV for the factor 
(MwiMP.nucieus)^ /^nucleus in © • Wc havc assumcd that 
the local WIMP density is 0.3 GeV/cm^, the WIMP ve- 
locity dispersion is 270 km/s and assumed that the es- 
cape velocity is infinite. By using that the Sun moves 



around the center of the galaxy with a velocity of 232 
km/s and that the Earth rotates around the Sun with a 
velocity of 30 km/s at 60 degrees inclination with respect 
to the galactic plane, we have found that the DAMA an- 
nual modulation amplitude implies (see ,18i. t24j for de- 
tails about such calculations): 
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Here ctwimp-simp is the WIMP-SIMP cross section, 
/simp is the abundance of anomalous iodine nuclei. If 
fwiMP-siMP is equal to Cdd in Eqn. ©, then: 
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This is less than the best limit of 6 x 10~^^ for light 
anomalous gold isotopes and far less than the limit of 
1 X 10"^ for heavy gold isotopes Of course, the 

lower the abundance of SIMPs, the higher ctwimp-simp 
can be before SIMPs would interfere with galactic struc- 
ture formation. So, the assumption ctwimp-simp = o^dd 
leading to |SJ| is very conservative. 

Most other direct detection experiments use elements 
that are lighter than iodine. The CDMS II experiment 
j25| uses germanium and the null results of this exper- 
iments pose the strongest constraints on spin indepen- 
dent WIMP-nucleon cross sections. The null results of 
GDMS II and the annual modulation observed by the 
DAMA experiment can be simultaneously explained if 
SIMPs don't bind to germanium or lighter elements and 
if the WIMP-nucleon interaction is sufficiently weak. The 
WIMP-nucleus cross section caused by the exchange of 
virtual SIMPs, cr^^n, should be of the order: 
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Here /j.wimp, nucleus is the relative mass of the WIMP- 
nucleus system, ^nucleus is the atomic mass number of 
the nucleus and g is the coupling between SIMPs and 
nucleons. According to the results of the CDMS II ex- 
periment, spin independent WIMP-nucleon cross sections 
larger than 4 x 10"'*'^ cm^ are excluded at 90% confi- 
dence level 25]. This corresponds to 2 x 10"'^^ cm^ for 
the WIMP-germanium cross section. We shall constrain 
the WIMP-SIMP model by demanding that cr™_„ < 1 
pico-barn for germanium. According to lO this implies: 
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Here /^wiMP,germ is the relative mass of the WIMP- 
germanium system and a; < 1 parametrizes the available 
window. Demanding that SIMPs don't bind to germa- 
nium nuclei constrains the coupling g. Using the assump- 
tion that SIMPs interact with nuclei in a similar way as 
hyperons, a relation between the minimum atomic mass 
number, Amin, below which SIMPs don't bind to nuclei. 
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FIG. 1: The upper and lower curves are Msimp and Mwimp 
as a function of logj^Q(l/x), respectively. 



and the SIMP-nucleon cross section was obtained in [15| . 
If we assume that the SIMP-nucleon cross section is 
barnes, this relation implies: 

9 ^2/3 ^ 1 , 1 (o\ 

2 GeV AfsiMP Af^„.„ ' ^ ' 

If we take ^min ~ 100, g becomes: 

. = 0.1(10- + -^). (9) 

Using the conditions Q, ^ and ©, A/wimp and 
.^siMP can be calculated as a function of the parame- 
ter X defined in Eq. [7| see Fig. ^ 

E.g., if a; = 10~^, a WIMP detector could detect nu- 
clear recoils in germanium consistent with a spin inde- 
pendent WIMP-nucleon cross sections of 2 x 10~^^ cm^ 
for a WIMP mass of about 50 GeV, apparently proving 
that the DAMA/Nal results are wrong. However, accord- 
ing to the WIMP-SIMP model and the assumptions we 
have made here, this would still be consistent with the 
DAM A results if A/wimp = 110 GeV and A/simp = 850 
GeV. 

Three other experiments are using elements heavier 
than iodine and have not detected nuclear recoils due to 
DM interactions. The CRESST II experiment 0| uses 
CaW04 and has a lower nuclear lower nuclear recoil de- 
tection threshold, 12 KeV for tungsten compared to 22 
KeV for iodine in the DAMA/Nal experiment. No recoil 
events in tungsten were detected during a 20.5 Kg day 
exposure. We have calculated that if the abundance of 
anomalous tungsten isotopes in the CRESST II detec- 
tor were the same as the abundance of anomalous iodine 
isotopes in the DAMA/Nal detector, there should have 
been about 17 events. This means that the abundance 
of anomalous tungsten isotopes should be at least an or- 
der of magnitude less than the abundance of anomalous 



iodine isotopes. Because SIMPs don't penetrate a large 
distance below the surface, the abundance of anomalous 
isotopes will strongly depend on the geological history of 
the material. Samples of tungsten and iodine typically 
have different geological histories |23|. Tungsten origi- 
nates from deep within the Earth. All tungsten deposits 
are of magmatic or hydro thermal origin. Iodine, on the 
other hand, is usually found in salt mines. The miner- 
als found in salt mines were originally dissolved in a sea. 
The salt mine formed when that sea receded and dried 
out. The SIMP-nucleon cross section imphed by © im- 
plies a stopping distance of the order of tens of kilome- 
ters So, unless the tungsten used in the CRESST 
II experiment came from very old deposits, the iodine 
in the DAMA/Nal detectors should have a larger abun- 
dance of anomalous isotopes than the tungsten used in 
the CRESST II experiment. Unfortunately, the CRESST 
II collaboration hasn't got the information about the lo- 
cation where their tungsten was obtained from 23>\. 

The DAMA/Xe-2 g9j and the ZEPLIN I [s^ exper- 
iments use liquid xenon. Since xenon has roughly the 
same atomic mass number as iodine, the same abundance 
of anomalously heavy isotopes would lead to the same ap- 
parent spin independent WIMP-nucleon cross section. In 
the DAMA/Xe-2 experiment xenon enriched in ^^^Xe was 
used. No events were detected and this was consistent 
with the results of the DAMA/Nal experiment assum- 
ing spin independent WIMP interactions. The ZEPLIN 
I experiment used natural xenon. The 90% exlusion 
limit from this experiment for spin independent WIMP- 
nucleon cross sections is about 1x10^^ pico-barns, about 
a factor of 2 below the lower limit of the DAMA/Nal ex- 
periment. Assuming that the WIMP-SIMP model is cor- 
rect, this means that the abundance of anomalous xenon 
isotopes in the ZEPLIN detector is less than the abun- 
dance of anomalous iodine isotopes in the DAMA/Nal 
detector. Like in case of tungsten, this could be due to 
a geological effect. There is some evidence that heavy 
noble gasses in the atmosphere are recycled in the man- 
tle It is thus possible that (very) heavy isotopes of 
xenon in the atmosphere are depleted by such a process. 

In the near future, the WIMP-SIMP model presented 
here can be verified or ruled out using data from dark 
matter direct detection experiments that are currently 
underway. If experiments using different materials report 
results that are consistent with a single WIMP model, 
then that would rule out the WIMP-SIMP model. If, on 
the other hand, only experiments using elements lighter 
than some critical atomic mass number are consistent 
with each other while experiments using heavier elements 
detect signals that are too large to be explained using 
conventional WIMP models, then that would be consis- 
tent with the WIMP-SIMP model. Experiments that 
use xenon are especially interesting. Using isotope sep- 
aration methods to enhance to concentration of lighter 
xenon isotopes one could deplete the concentration of 
anomalously heavy isotopes, thus influencing the WIMP 
detection rate. Since xenon has 9 stable isotopes of which 
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two have an odd atomic mass number this effect should 
be distinguishable from the predictions of any conven- 
tional WIMP model. 

In conclusion, we have shown that the annual mod- 
ulation signal detected by DAMA can be explained by 
WIMPs interacting with SIMPs captured in iodine nu- 
clei. The WIMP-SIMP cross section can be in the same 
range as is postulated for self-interacting dark matter 
without violating experimental bounds on anomalously 
heavy isotopes. 
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